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Star formation in molecular cores III. 
The effect of the turbulent power spectrum. 
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Abstract. We investigate the effect of the turbulent power spectrum (P(k) cc k^", with n = 3, 4 or 5) on the fragmentation 
of low-mass cores, by means of SPH simulations. We adopt initial density profiles and low levels of turbulence based on 
observation, and for each ii-value we conduct an ensemble of simulations with different initial seeds for the turbulent velocity 
field, so as to obtain reasonable statistics. We find that when power is concentrated at larger scales (i.e. for larger n), more 
protostellar objects form and there is a higher proportion of low-mass stars and brown dwarfs. This is in direct contrast with the 
recent results of Delgado Donate et al., presumably because they adopted much higher levels of turbulence. 
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1. Introduction 

Stars form in dense molecular cores (eg. Andre et al. 2000), 
and most stars - especially young stars - are in multiple sys- 
tems (e.g. Duquennoy & Mayor 1991; Mathieu 1994; Duchene 
1999, Patience et al. 2002), which implies that star-forming 
cores usually fragment into multiple objects. 

Simulations suggest that cores are prone to fragment into 
multiple objects under a variety of circumstances: (i) if they 
possess a small amount of initial rotation (eg. Burkert & 
Bodenheimer 1996); (ii) if their collapse is triggered by a sud- 
den increase in external pressure (eg. Hennebelle et al. 2003, 
2004); or (iii) if they contain turbulence. The level of turbu- 
lence can be high (eg. Bate et al. 2002, 2003; Delgado Donate 
et al. 2004) or low (eg. Goodwin et al. 2004a,b). 

This is the third in a series of papers investigating the 
collapse and fragmentation of cores with initial density pro- 
files and levels of turbulence based on observation. In Paper I 
(Goodwin et al. 2004a) we have shown that cores with even 
a very low level of turbulence can fragment into multiple ob- 
jects, but the number of fragments that form is very sensitive 
to the details of the initial turbulent velocity field. In paper 
II (Goodwin et al. 2004b) we have shown that the number of 
fragments increases as the level of turbulence is increased. In 
this paper we investigate the effect of the power spectrum of 
turbulence on the fragmentation of cores with low levels of 
turbulence. A similar investigation has already been made by 
Delgado Donate et al. (2004) for cores with high levels of tur- 
bulence, and there are significant differences between their re- 
sults and ours, which we explain in Section|4] 



In Section 2 we describe the initial conditions and nu- 
merical methods used. In Section 3 we present our results, in 
Section 4 we discuss them, and in Section 5 we summarise our 
main conclusions. 



2. Initial conditions and numerical method 

The density profiles of prestellar cores are approximately flat 
in the centre, and then decrease as r"^ with 2 < v < 5 in their 
outer parts, until they merge with the background (e.g. Ward- 
Thompson et al. 1994, 1999; Andre et al. 1996, 2000; Tafalla 
et al. 2004; Kirk et al. 2005). A good fit to the density profile is 
given by 
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where pkemei is the central density and ^kernel is the radius 
of the region in which the density is approximately uniform 
(cf. Whitworth & Ward-Thompson 2001). We set pkemei - 
3 X 10""* g cm""* and T^kernei - 5, 000 AU, with the outer bound- 
ary of the core at Rcok - 50, 000 AU, so the total mass of the 
core is M^ok = 5.4 Mq. These parameters are typical of low- 
mass star forming cores (eg. Jijina et al. 1999). The core is 
initially isothermal, with T = lOK, and molecular, with mean 
gas-particle mass fri - 4- x 10"^^ g. Hence the core has a ratio 
of thermal to gravitational energy of 
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The line widths of molecular cores show a significant non- 
thermal contribution (eg. Myers 1983; Myers et al. 1991; Jijina 
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Fig. 1. The filled circles give estimated values of atuib (the ratio of turbulent to gravitational energy) and Mcore (core mass) 
for the starless cores in the Jijina et al. (1999) catalogue. The open star shows the values used in this paper: atuib =0.1 and 
Mcore - 5.4Mq. The open circle shows the values used by Delgado Donate et al. (2004): cKturb = 1.0 and Mcore = S.OMq. The 
open square shows the values used by Bate et al. (2002, 2003): cKturb = 1 .0 and Mcore = 50Mq. Some of the cores in the Jijina et al. 
catalogue have akeady started to collapse, and therefore systematic infall motions are making a contribution to their non-thermal 
line-widths, thus causing us to over-estimate their o'tuib- Hence, we believe that our choice of ffturb is more representative of the 
initial conditions in these cores. 



et al. 1999), which is attributable to internal turbulence. Fig.[n 
shows the estimated ratios of turbulent to gravitational energy. 



and the estimated masses, Mcore, for prestellar cores from the 
Jijina et al. (1999) catalogue. These cores have been selected as 
prestellar on the basis of having low temperature (< 20 K), no 
associated IRAS source and no observed outflow. Using Fig.^ 
as a guide, we adopt ffturb = 0.1, and hence a total virial ratio 
of Cturb + fftherm ~ 0.4, for all the simulations in this paper. 

To model the turbulence in a core, we impose a divergence- 
free gaussian random velocity field with power spectrum 
P(k) oc r" (cf. Bate et al. 2002, 2003; Fisher 2004; Bonnell 
et al. 2003; Delgado-Donate et al. 2003, 2004). The observed 
velocity fields in GMCs and cores are well represented by tur- 
bulent power spectra of this form with « = 3 to 4 (Burkert & 
Bodenheimer 2000). 

In this paper we present two ensembles of 10 simulations 
each, one with n = 3 and one with « = 5, and compare these 
with the ensemble of 20 simulations with « = 4 already pre- 
sented in Paper II. Within an ensemble, each simulation differs 
only in the random seed used to initialise the turbulent velocity 
field. 

2.1. Computation method and constitutive physics 

The simulations are performed with the dragon code (Goodwin 
et al. 2004a), which is based on a standard implementation of 
SPH (eg. Monaghan 1992). An octal tree (Barnes & Hut 1984) 



is used to evaluate gravitational accelerations and to identify 
SPH neighbours. Particle smoothing lengths are adjusted so 
that each particle has A^neib = 50 + 5 neighbours. Gravity is 
kernel softened with the particle smoothing lengths, and stan- 
dard artificial viscosity is included with o-,. = 1 and yS,, = 2. The 
interested reader is referred to Paper I for further details. 

At low densities, radiative cooling is efficient and the gas 
in a core is approximately isothermal at r„ ~ 10 K. However, 
once the density exceeds pait ~ 10"'^ g cm"^, the optical depth 
through a core becomes too large for efficient cooling (Larson 
1969; TohHne 1982; Masunaga & Inutsuka 2000) and the gas 
switches to being approximately adiabatic. We model this be- 
haviour with a barotropic equation of state (eg. Tohline 1982; 
Masunaga & Inutsuka 2000): 
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Here P is the pressure, p is the density, and x 0.19kms ' is 
the isothermal sound speed in molecular gas at T ^ 10 K. 

Following the evolution of dense regions that are collaps- 
ing to stellar densities is very expensive computationally. In 
order to avoid this expense, wherever a bound regions forms 
with p > lOOpciit we replace it with a sink particle (Bate et al. 
1995). A sink particle interacts with the gas gravitationally, and 
it accretes any SPH particle which (a) approaches closer than 
10 AU and (b) is bound to it. When an SPH particle is accreted 
by a sink particle, the sink particle acquires the mass, linear 
momentum and angular momentum of the SPH particle, and 
therefore we can still monitor the conservation of these quan- 
tities as a check on the fidelity of the code. We refer to sink 
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particles generically as 'objects'; and then, more specifically, 
as 'stars' when the sink mass is greater than O.O8M0, and as 
'brown dwarfs' when the mass is lower than this. 

3. Results 

We have performed two ensembles of 10 simulations each, one 
ensemble using P{k) oc k^^ and one ensemble using P(k) oc k^^ . 
These are then compared with the ensemble of 20 simulations 
using P(k) oc k^'^ reported in Paper II. All simulations treat 
cores with o-tui-b - 0.10. A summary of the results is given in 
Table[2 The details of each simulation are presented in Table|2] 

3.1. The number of objects formed 

The fragmentation of a collapsing, mildly turbulent core usu- 
ally starts with the formation of a primary protostar surrounded 
by a rotating, but very unrelaxed, disc. The inflow of material 
onto the disc is very inhomogeneous and irregular, and as it 
joins the disc it causes spiral arms to develop. If these arms 
become sufficiently dense, they fragment to form secondary 
companions (e.g. Goodwin et al. 2004a,b; Gawryszczak et al. 
2006). If n is low, most of the turbulent energy is concentrated 
on small scales, but the resulting inhomogeneities in the inflow 
are of such low amplitude that the spiral perturbations they seed 
in the disc tend to be dissipated by the shear in the disc rather 
than being amplified by self-gravity. Conversely, if n is large, 
most of the turbulent energy is concentrated on large scales, 
and although the resulting inhomogeneities are again of low 
amplitude, they are of sufficiently large mass that the spiral per- 
turbations they seed in the disc have a better chance of being 
amplified by self-gravity and fragmenting into secondary com- 
panions. As n is increased from « = 3 to n = 5, more power 
is invested in large-scale turbulence, and therefore there tend to 
be more objects formed. When n = 3, A/obj = 3.7 + 1.4; but 
when n - 5, this increases to A/obj - 5.5 ± 3.0. 

We stress that the large variance on A/obj is because this is 
a chaotic process, and two simulations from the same ensem- 
ble (same aturb and n) can produce two vastly different sets of 
objects. For example one simulation from the ensemble with 
n = 4 produces a quadruple, a binary and four singles, whilst 
several others produce just one star (see Table|2l. Nonetheless, 
a generic trend is seen. 

3.2. Mass functions and companion probabilities 

Once an object forms, its final mass is determined by compet- 
itive accretion (Bonnell et al. 2001) and dynamical interaction 
with other objects. Competitive accretion causes the more mas- 
sive objects, and/or those which reside in the dense material at 
the centre of the core, to grow rapidly in mass. Dynamical in- 
teraction causes some objects, usually the lower-mass ones, to 
be ejected from the dense material at the centre of the core, 
so that the remaining objects become more tightly bound and 
some eventually end up in stable multiple systems. Thus, in 
general, it is the lower-mass objects (low-mass stars and brown 
dwarfs) which are ejected as singles before they can accrete 



much mass (Reipurth & Clarke 2001); and the higher-mass ob- 
jects (~ 1 M0) which remain in the centre of the core, and form 
multiples. 

Figure |2] shows the mass functions from the three ensem- 
bles with n - 3 (top), n - A (middle) and n - 5 (bottom). In 
each case the filled portion of the histogram represents objects 
in multiple systems, and the open portion represents single ob- 
jects. 

The mass functions for n - A and n - 5 are very similar 
There is a broad peak around 1 M©, consisting mainly of ob- 
jects in multiple systems, and a flat tail of lower-mass objects 
consisting mainly of ejected singles. Because of the large num- 
ber of ejected singles, the overall companion probability is low, 
-0.6. 

In contrast, when n = 3, the mass function is dominated by 
the peak around ~ 1 and there are very few ejected singles. 
This is because, when « = 3, an individual core spawns fewer 
objects, and therefore fewer ejections are required to stabilise 
the remaining multiple system. The paucity of ejected singles 
gives a much higher overall companion probability, ~ 0.8. 

3.3. Binary separations and mass ratios 

In all cases (n = 3, 4 and 5), the distribution of separations is 
much narrower than that observed by Duquennoy & Mayor for 
local G dwarfs, which is not surprising, since we have consid- 
ered only one core mass, and only one level of turbulence'. In 
particular there is a total lack of wide binaries (a > 100 AU). 
For « = 3, most binaries have separations in the range 10 
to 30 AU, and there is only one hard binary (a < 10 AU). 
In contrast, when n - 4, 17 of the 20 binaries formed have 
a < 10 AU; and when « = 5, 10 of the 13 binaries formed have 
a < 10 AU. This difference arises because a core with larger n 
tends to spawn more objects, and so on average more ejections 
are needed before a stable multiple is left; specifically, the av- 
erage number of objects ejected from a core is only 0.80 for 
n - 3, but 1.95 for n - A and 2.10 for n - 5. Since each ejec- 
tion hardens the multiple that is left behind, the greater number 
of ejections for n - A and n — 5 means harder multiples, i.e. 
smaller separations. 

Fig. |3] shows the distribution of mass ratio, q = MJM^, 
against separation, a (strictly, semi-major axis). In all cases 
the mean mass ratio is high, viz. q - 0.1 A, 0.83 and 0.65 for 
n = 3, 4 and 5, respectively. There is a tendency for closer bi- 
naries to have higher mass-ratios, i.e. more nearly equal com- 
ponents. This tendency arises because the formation of a close 
binary often entails hardening by ejection, and the ejections 
tend to remove the less massive objects; hence only the more 
massive objects are left as potential binary components, and the 
range of possible masses is thereby reduced, pushing q towards 
unity. In addition, when a close binary accretes material with 
high angular momentum, the accreted material tends to end up 
on the less massive component, which again pushes q towards 



' Hubber & Whitworth (2005) have shown how the full range of 
binary parameters can be reproduced by considering a range of core 
parameters 
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Fig. 2. Normalised mass functions for objects spawned by cores with n - 3 (top), n = 4 (middle) and n - 5 (bottom). The filled 
portion of each histogram represents objects in multiple systems, while the open portion represents single objects. 
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Table 1. For each value of n, we list the number of realisations simulated (A/'ieai); the mean mass in objects at the end of a 
simulation (M^oi/Mq); the mean number of objects formed per simulation (^Nobj) and its variance; the mean number of stars 
ejected from a core (A/'ej); the net ratio of brown dwarfs to stars (A/bd/A/^,); the fraction of all objects that are single (S /Nobj)', 
the median semi-major axis (flmed); the mean mass-ratio for binary systems (q); the variance of the mass ratio of binary systems 
{(Tq); and the mean numbers of singles (S), binaries (B), triples (T), quadruples (Q), and quintuples (Q') formed by one core. 
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Fig. 3. Mass ratio, q, against separation, a, for binaries spawned by cores with n - 3 (crosses), n - A (open circles) and n - 5 
(solid stars). 



unity (Whitworth et al. 1995, Bate & Bonnell 1997, Papers I & 
II). 

Such high mass ratios are not compatible with the observa- 
tions, which, whilst they show a trend to more equal-mass com- 
panions in close binaries (e.g. Mazeh et al. 1992; White & Ghez 
2001; Fisher et al. 2005), are certainly not as extreme as sug- 
gested by these results. However, a recent paper by Yasuhiro 
et al. (2005) suggests that accretion in a proto-binary system 
is generally onto the primary (especially in circular orbits) as 
angular momentum is removed from the accreting gas by spi- 
ral shocks. Such a mechanism is beyond the ability of these 
simulations to resolve and may help solve this problem. 

4. Discussion 

Delgardo Donate et al. (2004) have also explored the effect of 
the turbulent power spectrum on the fragmentation of low-mass 
cores, performing two ensembles of 5 SPH simulations each, 
one with n = 3 and one with n - 5. They use a slightly dif- 
ferent core mass (5.0 Mq, as compared with our 5.4 Mq), a dif- 
ferent density profile (uniform, as compared with Eqn.Q, and 
a slightly different equation of state. However, the most signif- 



icant difference is that their cores have a much higher initial 
level of turbulence; specifically, they adopt Cturb = 1-0, as com- 
pared with our aturb = 0.1. As a consequence, the initial turbu- 
lent velocities in their cores are mildly supersonic and the cores 
are marginally unbound (aturb + Q'theim - 1-1), whereas the ini- 
tial turbulent velocities in our cores are subsonic and the cores 
are approximately virialised (o-turb + Q'therm - 0.4). The amount 
of turbulent energy has a profound influence on the outcome of 
collapse. 

In the strongly turbulent cores of Delgado Donate et al. 
(2004), there is so much power in the turbulence that even the 
small-scale inhomogeneities created by small-scale turbulent 
motions can become self-gravitating and collapse (whereas in 
our simulations these small-scale inhomogeneities have much 
lower amplitude and tend to disperse). As a result, many more 
objects are formed, and - in direct contrast with our results 
- more objects are formed when n = 3 than when n - 5. 
This is because the same amount of power invested in small 
scales (large k) produces a larger number of inhomogeneities 
than when it is invested in large scales (small k). Furthermore, 
if more objects are formed in a core, then there have to be more 
ejections before a stable multiple is created, so the mass func- 
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Table 2. The results of the individual simulations; in all cases fftuit - 0.10. Column 1 gives the simulation identifier and Column 
2 gives n (the exponent of the turbulent power spectrum, P(k) oc k^"). Column 3 gives Mobj (the total mass of objects formed, 
stars plus brown dwarfs). Column 4 gives A/'obj (the total number of objects formed) and Column 5 gives A/bd (the total number 
of brown dwarfs formed). Column 6 gives the multiplicities of the multiple systems formed, and Column 7 gives the mass of 
each individual object. Those objects which are part of a binary system are distinguished with those which are part of a triple 
system with ', and those which are part of a quadruple or quintuple system with Two realisations have ejected binary systems. 



ID 


n 






A'bd 


Multiplicity 


Masses/Mo 


571 


3 


3.04 


4 





Triple 


1.06', 0.77', 0.61,0.59' 


572 


3 


3.21 


5 





Quadruple 


1.24'', 0.75^ 0.60', 0.47, O.W 


573 


3 


3.62 


5 


1 


Quadruple 


1.43*, 0.78', 0.75', 0.62', 0.04 


574 


3 


3.20 


5 





Triple 


0.89', 0.80', 0.80', 0.57, 0.14 


575 


3 


2.95 


1 





Single 


2.95 


576 


3 


3.79 


5 


1 


Quadruple 


1.05', 1.05', 0.86', 0.74', 0.09 


577 


3 


3.58 


2 





Binary 


2.01*, 1.57* 


578 


3 


2.97 


4 


1 


Triple 


2.04', 0.45', 0.45', 0.02 


579 


3 


3.79 


2 





Binary 


2.47*, 1.33* 


580 


3 


3.76 


4 





Quadruple 


1.36', 0.81', 0.80', 0.79' 


001 


4 


3.78 


3 





Triple 


1.49', 1.15', 1.13' 


002 


4 


2.83 


1 





Single 


2.38 


003 


4 


3.72 


1 





Single 


3.72 


004 


4 


3.48 


1 





Single 


3.48 


005 


4 


2.86 


4 


1 


Binary 


1.43*, 0.77, 0.65*, 0.02 


006 


4 


2.84 


1 





Single 


2.84 


007 


4 


3.15 


5 





Triple & Binary 


1.76', 0.72', 0.47', 0.10*, 0.10* 


008 


4 


3.22 


6 


2 


Quadruple 


1.97', 0.47', 0.35', 0.34', 0.06, 0.03 


009 


4 


3.48 


8 


4 


Quadruple 


2.28', 0.49', 0.26', 0.25', 0.08, 0.05, 0.04, 0.04 


010 


4 


3.31 


8 


1 


Quadruple 


0.76', 0.74', 0.58'', 0.57', 0.46, 0.09, 0.08, 0.03 


Oil 


4 


3.96 


12 


4 


Triple & binary? 


0.89', 0.82', 0.82', 0.42, 0.38, 0.25, 0.12, 0.11, 0.04*, 0.04*, 0.03, 0.03 


012 


4 


3.60 


6 


2 


Triple 


1.34', 0.92', 0.79', 0.50, 0.04, 0.02 


013 


4 


3.18 


10 


3 


Quadruple & binary 


0.77', 0.68', 0.61', 0.60«, 0.11*, 0.11*, 0.10, 0.06, 0.05, 0.04 


014 


4 


3.29 


4 


1 


Binary 


1.58*, 1.16*, 0.49, 0.08 


015 


4 


2.48 


1 





Single 


2.48 


016 


4 


3.58 


4 





Triple 


1.23', 1.15', 1.11', 0.09 


017 


4 


3.41 


8 





Quintuple 


1.10', 0.98', 0.32', 0.27«, 0.27', 0.17, 0.15, 0.14 


018 


4 


3.48 


4 





Quadruple 


0.98', 0.94', 0.79', 0.77' 


019 


4 


3.58 


5 


1 


Triple 


1.38, 1.03, 1.00,0.11,0.06 


020 


4 


3.77 


3 





Triple 


1.28', 1.27', 1.22' 


551 


5 


3.21 


1 





Single 


3.21 


552 


5 


3.34 


7 





Quadruple 


0.69", 0.59', 0.57', 0.52", 0.45, 0.37, 0.15 


553 


5 


3.46 


6 





Triple 


1.29', 1.11', 0.49, 0.29', 0.17, 0.10 


554 


5 


3.22 


8 


2 


Quadruple 


1.10', 0.75', 0.71', 0.27', 0.18, 0.15, 0.04, 0.01 


555 


5 


3.52 


5 





Quintuple 


1.16', 0.84', 0.60", 0.51', 0.40' 


556 


5 


3.39 


1 





Single 


3.39 


557 


5 


3.77 


7 


1 


Triple + Binary 


1.45', 0.93', 0.89', 0.22*, 0.13*, 0.13, 0.01 


558 


5 


3.10 


11 


3 


Quintuple + Binary 


0.83", 0.57', 0.43', 0.35', 0.30', 0.29, 0.14, 0.09, 0.04, 0.02*, 0.02* 


559 


5 


3.75 


6 


2 


Triple 


1.13', 0.97', 0.80', 0.79, 0.03, 0.03 


560 


5 


3.30 


3 





Triple 


2.25', 0.59', 0.46' 



tions derived by Delgado Donate et al. (2004) have a larger tail thermal linewidth is attributable to collapse rather than initial 
of low-mass singles (low-mass stars and brown dwarfs), and turbulence. The low levels of turbulence we have invoked also 
this tail is larger when n - 3 than when n - 5. seem to result in more acceptable values for the overall multi- 

With reference to Fig. ^ we suggest that the level of tur- pHcity. 
bulence we have adopted is more representative of the cores in 
the Jijina et al. (1999) catalogue, and therefore probably more 
representative of the initial conditions in low-mass star forming 

cores, particularly if, as seems likely, some of the Jijina et al. We have investigated the influence of the slope of the turbulent 



5. Conclusions 

We have investigate 

cores are already collapsing, and therefore some of their non- power spectrum on the fragmentation of dense molecular cores, 
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by means of an ensemble of SPH simulations. We consider 
a spherical, 5.4Mq core, with a Plummer-like density profile 
(Eqn.0, and a low level of turbulence, atuib = £^tuib 

/|Q| = 0.10, 

similar to observed cores such as L1544. The turbulence has a 
power spectrum P(k) oc k^" with n = 3, 4, or 5. The choice of 
n influences the number of objects that form, the mass function 
of those objects, and the properties of the multiple systems that 
they comprise. However, the process is chaotic, in the sense 
that even if n is fixed, different realizations of the turbulent ve- 
locity field can produce widely different stellar masses and bi- 
nary properties. Our main conclusions are therefore statistical 
in nature, and could always be improved by performing more 
simulations: 

- The average number of objects that form in a collapsing 
core increases monotonically with n, from A/'obj = 3.7 + 1 .4 
when « = 3, to A/'obj = 5.5 + 3.0 when n - 5. 

- The mass function always involves a peak at M ~ IM© and 
most of the objects in the peak are in multiple systems. 

- As n increases and more objects are produced, more dy- 
namical ejections are required before a stable multiple is 
formed, which has two consequences, (i) A larger propor- 
tion of single, low-mass objects (low-mass stars and brown 
dwarfs) is produced; hence the mass function develops a 
low-mass tail and the mean multiplicity decreases, (ii) The 
resulting binaries tend to be harder, i.e. to have smaller sep- 
arations. 

- The mean mass ratios of binary systems do not depend 
strongly on «, but close binaries tend to have mass ratios 
closer to unity, i.e. more nearly equal components. This is 
because the ejections which harden a binary preferentially 
remove low-mass objects, leaving the two most massive ob- 
jects; the larger the initial number of objects, the more ejec- 
tions are required, the harder the final binary, and the closer 
the masses of the two binary components. 

- The low level of turbulence we have adopted in these sim- 
ulations (atuib =0.1) appears to be in good agreement with 
observations of low-mass cores, and to reproduce the mean 
multiplicity in observed stellar populations. 

These results do not match observed multiple systems very 
well. The separation distribution is too narrow and the mass 
ratios tend too much towards equal-masses. This is probably 
due in part to the fact that we consider only one (rather high) 
core mass. The separation distribution may well improve if a 
realistic range of core masses were considered (cf. Hubber & 
Whitworth 2005), but this would vastly increase the compu- 
tational load as a far larger series of ensembles would be re- 
quired. Even then, the small separation range is exacerbated 
by the hardening of systems through ejections when too many 
fragments are formed (see Goodwin & Kroupa 2005). The ten- 
dency to equal-mass binaries may be a result of unresolved 
physics in the inner accretion region (see Yasuhiro et al. 2005). 
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